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1 6 Abstract 


A preliminary study to compute the size of a special 
radiometric reflector antenna is presented. Operating at 
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1. Sununary 


A preliminary study to compute the size of a special 
radiometric reflector antenna is presented. Operating at 
1 GHz, this reflector is required to produce 200 simul- 
taneous contiguous beams, each with a 3 dB footprint of 
1 km from an assumed satellite height of 650 km. The 
overall beam efficiency for each beam is required to be 
more than 90%. 
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s 2. Introduction 

\ 


The purpose of this report is to present a pre- 
liminary study for tht design of a large radiometric 
reflector antenna system. When orbiting at a height 
of 650 km, this antenna system is required to produce 
simultaneously 200 contiguous 3 dB circular footprints 
on the ground, each having a diameter of 1 km. The 
lowest frequency of operation is 1 GHz. The footprints 
are required to be as identical to each other as possible. 
The single most important requirement on the system is 
that the overall beam efficiency for the copolarized 
component in each of the 200 beams be better than 90% 
within the two and a half 3 dB beamwidths. This means 
that among other things, the cross polarization be 
minimum too (<25 dB) . 

When in orbit, the reflector may undergo considerable 
thermal distortions and its performance may change. A 
method is therefore needed to predict the performance of 
even a distorted reflector. Such a technique is discussed 
in Section 6 of this report. 


a.' 




3* Spherical Reflector Approach 


A solution that meets the requirements set in Section 
2 is schematically shown in Figure 1. The 200 beams are 
simultaneously obtained by stacking 200 identical feed 
ant3nnas along a concentric circular arc in front of a 
spherical reflector such that each feed is pointing 
radially towards the spherical reflector surface. Each 
feed thus creates its own independent footprint. And 
since each feed antenna essentially sees an identical 
segment of the spherical reflector, the resulting 200 
footprints are also practically identical. Observe that 
the angular separation between any two consecutive feed 
antennas (called B) is the same as the angular separation 
between the two adjacent footprints. This, for an altitude 
of 650 km and a footprint size of 1 km, turns out to be 
0.088® and the 200 feeds stacked along the feed arc thus 
subtend a total of 16.6° angle at the center of the spherical 
reflector (Figure 2). 

Notice that the angular separation between any two 
consecutive feeds depends only upon the altitude and the 
footprint size; the physical separtion, however, is the 
product of the angular separation (0 = 0.088°) and the 
feed arc radius, and therefore, depends upon the radius 
of the feed arc also. The radius of the feed arc, there- 
fore, should be large enough to provide enough physical 
room for each of the 200 feed antennas. It is assumed 
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Figure 1 — Multibeam Spherical Reflector Concept 
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Figure 2 — Geometry of Sphe, }cal Reflector Antenna 


ft thi# stage t)^^t each €jse4 antenna should have# on |n 
average, a rgcp of et least 88 cm, which leads to e 
arc radios of 577 m. And since for spherical ref^ctoirpif 
the fjsed arc is generally located about halfway bajb^eeiii 
the reflector end its center of curvature (nearer to the 
reflector) , the radigs of curvature of the spherical 
reflector is chosen to be 1175 m. Note that there is go 
specific reason tQ pick 88 cm for feed spacing except fhat 
the estimations of practical feed horn sizes suggest thg^ 
a room of about a meter bo available for each feed antenna* 

And of course, whether a feed horn limited in size to 88 cm 
at 1 GHz feeding a reflector with dimensions chosen above 
can give a satisfactory secondary far field pattern or not, 
remains to be checked. 

Let us now consider an individual footprint which is 
caused by an individual feed antenna located at the feed 
arc. Each feed antenna illuminates a portion of the spherical 
reflector and it is the far field of this illuminated reflec- 
tor aperture which must have (a) a 3 dB beamwidth of 0.088**, 
and (b) a beam efficiency of better than 90% within the two 
and a half 3 dB beamwidths. The later requires that the 
highest side lobe of the reflected pattern be less than -32 
dB with wide angle side lobes below -80 dB. A study of 
various aperture distributions [1] indicates that for an 
operating frequency of 1 GHz, an aperture d tameter of about 
300 m (say, 305 m) with a rotationally symmetric cosine 
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squared field distribution produoee both a 3 dB beanwidth 
of 0.088* and a side lobe at -32 dB# the side lobe fall off 
being -18 dB/Octave* 

Returning to Figure 2# an illuminated aperture 
with a diameter of 305 m on the reflector corresponds to a 
cone of 15* half angle emanating from each feed antenna. 
Therefore# each feed antenna whose nominal diameter has been 
fixed to 88 cm at 1 GHz has to be able to produce a rotation- 
ally symmetric cosine squared far field pattern over + 15*. 
The overall diameter of the spherical reflector dish to 
produce 200 beams then turns out to be 660 m as shovm in 
Figure 2. 


• • 


♦ • F— d Coni Iterations 

In the previous section# it was asstoned that after 
reflection# each feed pattern gave rise to a rotationally 
symmetric aperture field which varied as cosine squared in 
the radial direction. For the reflector dimensions unde.: 
consideration# the portion of the reflector illuminated by 
a feed is such a small fraction of the full sphere# that 
it is practically flat and therefore a rotationally symmetric 
cosine squared aperture distribution is easily achieved by 
a feed which too has a rotationally synnnetric cosine squared 
pattern. 

A rotationally Sj ranetric cosine squared feed pattern 
can be generated by any one of the several types of horns. 

In the present study# however, a circular corrugated horn [21 
is considered. The feed pattern of the corrugated horn used 
for the computations presented in the following sections is 
shown in Figure 3. Observe that the horn diameter at the 
mouth is 2 m which is larger than 88 cm, the space designated 
for each feed at the feed arc. Therefore# the '.eed horns 
will have to be staggered around the feed arc so that they 
are still on an average 88 cm or 0.088® apart. The beam 
efficiency of the feed pattern within 15", which corresponds 
to the edge of the reflector is 98.3%. This number is 
imporcant because the overall beam efficiency of the antenna 
system is the product of this efficiency and the beam efficiency 
of the secondary pattern. 
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Figure 3 




s. 


Secondary Pattern 


For the ref Ic geometry shown in Figure 2 and 

using the feed pattern presented in Section 4, the computed [3] 
secondary radiation pattern is shown in Figure 4. It haa 
a 3 dB beiitiiwldth of 0.080** and a maximum cross polarisation 
level of less than -200 dB. The beam efficiency of the 
secondary pattern at two and a naif 3 dB beamwidths is 93.4%, 
the overall beam efficiency, therefore, being better than 91%. 

one of the concerns in spherical reflector applications 
is the resulting spherical aberration, it is of interest to 
note that in the present case, such a small segment of the 
sphere is being used as reflector that the maximtim spherical 
aberration near the edge of the illuminated aperture (where 
the field strength is -26 7 dB, Figure 3) is equivalent to 
a phase error of only about 18®. Such a small phase error 
causes a negligible degradation in the antenna gain. 


original is 
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Angie off the Axis in Degrees 


6. Thermal Distortion Considerations 

The performance of a reflector antenna in space is 
sometimes not the same as predicted by the initial design 
because the reflector undergoes severe distortions due to 
thermal variations. If the distorted shape of the reflector 
is quite a bit different fr<wi the original spherical shape, 
the reflector performance may change significantly and may 
even become unsatisfactory. Therefore, it will be desirable 
to be able to predict the performance of even the distorted 
reflector. If the distorted reflector surface could be knoWn 
analytically, then the reflector performance ofcourse could 
be accurately predicted. It is not generally possible to 
know an analytic expression for the entire distorted reflector 
surface at all times. Alternatively, a sampling scheme can 
be implemented such that the coordinates of many discrete 
target points located along a rectangular grid on the reflector 
surface are known. Then, a smooth tight cubic surface can 
be fitted through the four corners of each of the rectangular 
grid patches such that the whole composite reflector surface 
is continuous and has continuous partial derivatives. Using 
this piecewise analytic expression for the reflector surface, 
the reflector properties can be computed. Needless to say, 
the target points on the reflector surface must be dense 
enough to sample the distortions and such that the surface 
between the measured points could be assumed tightly stretched. 



To demonstrate that the far field radiation pattern can 
indeed be accurately computed even when the reflector surface 
is known only at certain discrete points, the following 
example is presented. 

Computations presented in Figure 4 are made again, 
this time, instead of using a single analytic expression for 
the entire spherical reflector surface, though, the x-, y-, 

2 “ coordinates of 45 equispaced points on the reflector 
surface are used. These surface points are located on the 
reflector surface along a rectangular grid as shown in 
Figure 5, the points being 40 m or 133.3 wavelengths apart. 

For computational purposes, the reflector surface over 
any rectangular patch is expressed as a bi-spline under 
tension [4]. In Figure 6, the far field radiation pattern 
computed by using a single analytic expression for the entire 
spherical reflector surface (as in Figure 4) is shown by 
solid lines. On the same figure, the far field radiation 
pattern computed by using the piecewise analytic composite 
surface through 45 target points on the reflector is plotted 
with solid dots. The field values not shown by solid dots 
are too close to the solid line curve values to distinguish. 
The modifications needed in the computer program in Reference 
3 to make the present surface fitting computations are shown 
in the Appendix, 

In conclusion, for actual distorted reflector conditions, 
where the whole distorted reflector surface is not known 
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Figure 5 — Location of Surface Target toints used for Surface FiWiiig Catcufations 
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Figure 6 — Computed Secondary Pattern 


analytically » accurate far field c<Hiiputatlona can be maAn 
by uaing our computer program which will accept for the 
reflector surface geometry » a set of discrete reflector 
surface points. The basic underlying asstmiption is that 
;.he surface is smooth between the target san^le points. 
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Appendix 

The computer program documented in Reference 3 haa been 
modified auoh that now It can be used for making alao the aurfaof 
fitting reflector calculations similar to the ones presented in 
Section The purpose of this appendix Is to document the oerres* 
ponding modifications* changes* and additions. 

T%ro new features have been added to the computer program msFX*CTR 
{documented In Reference 3). The first feature* which is not of 
direct concern to the subject matter of this report is that In 
addition to par 2 d>olic, spherical* and ellipsoidal reflectors* the 
program REFLCTR can now also handle planar reflectors. A planar 
reflector is specified by (a) three cartesian coordinates of a 
point on the surface of the planar reflector - PLNPNT(l)* PLNPNT(2)* 
and PLNPNT(3), and (b) the three cartesian components of a unit 
vector normal to the reflector surface - PLNORM(l), PLN0RM(2)* and 
PLN0RM(3). The value of the integer variable SURFACE must be set 
to 4 for a planar reflector. 

The second new feature is that in addition to specifying a 
reflector surface as being parabolic, spherical, ellipsoidal, or 
planar by setting SURFACE = 1,2,3 or 4 respectively, one can now 
also choose to specify any of the above reflector surfaces in 
terms of only a finite number of discrete target points located 
along a rectangular grid (assumed square grid here for simplicity) 
on the reflector surface. This is done by setting the integer 
variable NEEDFIT to a nonzero positive value. SIGMA is the tension 
factor (defined later) used for fitting the surface through the 
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grid points f the grid spacing being DISFAC for both the y<- 
and the z~ directions. All these newly defined variables 
along with the ones already defined in Reference 3 are read# 
as before, in the subroutine NPUT. A listing of modified 
NPUT is given in Figure A-1. 

As a result of implementing the above two features, the 
subroutine APERTUR also changes. A listing of the new APERTUR 
is given in Figure A-2. Notice that in this subroutine, 
before statement number 100, the coordinates of surface 
target points are first stored in dimensioned arrays called 
EXTRAX, EXTRAY, and EXTRAZ and then the subroutine SURFl is 
called to compute the parameters necessary to compute an 
interpolatory surface passing through the surface grid 
points. Later on in the subroutine APERTUR (before statement 
number 130), subroutine SURFD2 is called to interpolate 
the reflector surface at the given coordinate pair and to 
compute the components of a normal vector at the inter- 
polated point. 

The subroutines SURFl and SURFD2 are from "A Spline 
Under Tension Package for Curve and Surface Fitting" by 
A. K. Cline, Department of Computer Science, University of 
Texas, Austin. This package of subroutines is an extension 
of Cline's work reported in Reference 4. A listing of 
subroutines SURFl and SURFD2 which also includes definition 
of parameters used in the subroutines is presented in 
Figure A-3. 



no 

140 

no 

160 

170 


200 

555 

556 
S7« 


57<» 


SIIRi^nuTlNF »<PI»T(P» 

CriMMOM/PAIAMS/TITl.FCn) ,AnRnRF,XLAM.OmD.SH«FACF.*PPnTA,FFFnm . 
, AI.PHA.OFTA.OAmma.XC ,YC.2C ♦HFMAcx^HFMIFX.RMTP.RMPP, 

. NT.NP,NPniMT.MAXPTS,RFl.l.P.»LMPWT(3) .PLNnRM(^) 

CnMMnM/PLnCKG/YCHI..7CHL.HFMA«l.,HHM!H(, • 
r.nMMON/PATTFRN/PHl ( ^> ,TH(-TA( ? ) 
r.nMMOM/MATH/Pt ,Pi2,Pio2,OTn^.RTnn 
COMMON/SUH F I T/NFFOF I T , S ! OMA , 01 SF AC 
IMTFGFR TI TLF.SlIRFACF 
BFAO 200, TtTt.F 

RFAO 4 . AnRnRF,XLAM,G‘Un,SlUFACF,MFFnFIT,APROTA 


RFAO 4, RFU.P .P'.MPNTfPLNORM.SlGMA ,niSFAC 
RFAO 4, FFFO, ALPHA, OFTA, gamma 

RFAO 4, XC,YC,ZC,HFMAFX,HFMIFX,YCRL,7CPL.HFMARL,HFMIBI. 

RFAO 4, PHI.THFTA 

IF( APROTA.GT.0.0) WR I TF ( ?0 , SRR ) T 1 Tl F 

IF( APROTA.GT, 0.0» WR I TF( 20, «5S6 ) pcpo, AOROR F , XLAM ,GR f 0, ALPHA , RFTA , TAPF *0 

. . GAM,.*A,X-:,YC,7C,HFMAFX,HFMrFX TAPF 20 


GO TO ( no, 140, no, >60) ,SMRF4^F 

PRINT 57B, TITLF,XLAM,FFF0, alpha, RFTA, GAMMA, PLNPMT.PLNORM 
GO TO 170 
PRINT 579, 

GO TO 170 
PRINT 5R0, 

GO TO 170 
PRINT 5«1, 

PRINT 5«2, 


4******* 


TI TLF,XL am, fffO, ALPHA ,HFT 4, gamma ,A nRORF,BFLLP 

TITLF,XLAM,FFF0,4LPhA,RFTA ,GAMMA ,A nRORF 

TITI.F , XL am, FFFO, alpha ,«FT A, gamma .AORORF 
XC,YC,7.C,HFMAFX,HFmIFX,GRI0, YCRL,7CRL,HFmARL,HFMIRL, 
THFTA,PHI 

PRINT 5RR. sigma, orSFAC 


IF(NFFOFIT.GT.O) 
formatiraio) 

FORMAT! 1X,RA10) 

FnRMAT(lX,F15.4) 

FORMAT! IHI ,////, nX. 4PLAMAR RFFI.FCTOR FAR FIFI.O RAOIATION •/ 

4 PATTFRNi CAI.CIILAT IPN4///4 4RA10/4 4RA1 

0/4 WAVFI.FNGTH OF FLFCTOC FIFI.O 4F9.4/ 

4 LOCATIOm of COOROINATF ORIGIN wR T FFFn (X.Y,7) 43F7.2 

/4 FFFO ROTATION AnglfS ! Ai.PHA , RctA .GAMma ) 4AF7.2 

/4 A POINT ON THF rffLFCTHR SlIRFACF (X,Y,7) 4RF7.2 

/4 COMPONFNTS OF UNIT Npi^HAl. TO SliRFACF (X,Y,7)-. 4RF7.2 


) 

FORMAT! IHl ,////. I IX. 4FI.LIPTir.AI. RFFI.FCTOR FAR FIFI.O RAOIATION */ 
4 RATTFRM calculation*///* 4RA10/4 4RA1 

0/4 WAVFI.FNGTH OF FLFCT-UC FIFI.O 4F9.4/ 

4 location of COOROINATF TRIGIN wRT FFFO |X,Y,7) 4RF7.2 

/* FFFO ROTATION ANGLFS { A I. PH A , R F T A , G AMM A ) 4RF7..? 

/4 MAJOR AXIS OF THF FI.L^TICAL RFFLFCTOR .4F7.7/ 


Figure A-1 
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5<»0 




5B? 


5<50 


600 


. * MINOB axis op THF PIUPTICAI. RFFLFCTnu. ...*F7.2) 

FORMAT! IHl ,////, llX.OSPMpqiCAI. RFFI.CCTOR FAR FIF(.0 RAOIATinN */ 
« PATTcRM calculation*///* *AA10/* *«A1 

n/* wavflfmgth of flfct^ic fifi.o ,.*fb,*/ 

* LOCATIOm of CnORDlwATF ORIGIN uRT FFFD ( X .Y * Z » , . . . ,*BF7.2 

/* FFFO ROTATinw AMGI.F' ( ALPHA , RFT A , GAMMA ) *BF7.2 

/* RAOIMS OF THF RCFLFCTOR SPHFRF. ,.,..*F7.2) 

FORMATIlHl,////, 11X4*PARAH0LIC RFFI.FCTOR FAR FIFI.O RAOIATION •/ 

* PATTFHN calculation*///* *«A10/* *RAl 

0/* WAVFI.FNGTH OF FI.FCTUC FIFI.0 ..*F«».«/ 

* LOCATIOM ilF COOROIMATP ORIGIN uR T FFFO IX. Y, 7) *BF7.2 

/* FFFO ROTATION ANGLES , A I.PHA , RFTA , GAMMA ) *:»F7.2 

/* FOCAL LFNGTH of THF RFPI.FCTOP *F7.2) 

FORMAT(=:- APFRTIIRF PI.ANF LOCATIOM ( xC ) * FR.2/ 

* CnOROINATFS OF THF APPRTHRF PI.ANF CFNTFR *2F7.2 

/* HALF MAJOR AXIS 0F‘ A'PCRTMRF pi.ANF (ALONG V *F7, 2/ 

* HALF minor axis of APFRTIIRF PI.ANF (ALONG 7 .*F7,2/ 

• » GRin SI7F ilSFO FOR mmmfRTCAI. INTFGRATION *F9.A/ 

* FFFO SHAOnw CFNTFR fnnRniNATFS IN APFRTIIRF PL *2F7,2 

/* half major axis of FFFO SHAOOH *F7.2/ 

* half minor axis of FFFO SHAOOW *F7,2/ 

»■ ThFTA RANGF FOR FFFO PATTFRN (L.H.I - OFGRFFS) *BF7,2 

/* PHI RAMGF for FFFO PATTFRN (L.H.I - OFGRFFS), *BF7,2 

) 

FORMAT!* SIGMA IlSFO FOR SliRFACF FITTING,,.., *F7,2/ 

. * SPACING BFTWFFN SiUFACF POINTS *F7,2) 


FORMAT!* .IM-SiiFFICIFNT WORK STORAGF, NFFOFO *IS* AwAILABLF IS 

ONLY *15* *) 

PI=AC0S(-1 .0) 

PI?=PI+PI 

Pin?=:0.5*PI 

OTOR=PI/1RO. 

RTnn=l*0./PI 

NP = ( PHI ( ?)-PHI { I ) I/i^Hl ( X) + l,5 

NT=( THFTA(7)-THFTA( 1) )/THFTA( B) + l .5 

NPOINT=NT*NP 

IFINPOINT .GF. maXPTS) go to 600 
CALL FILLIP. NT, NP) 

RFTtiRN 

PRINT 590, NPOINT. MAXPTS 

STOP 

END 


Figure A-1 (Continued) 
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SURUnilTINF APF«TIU( P.NTX.NPX) 

CnHMnM/PA«AMS/TITLE( 1 A U AORtU': , XI. A «,r,R I SMR FACP t APPOf A , PPPOI 1 1 « . 

. At.PHA,KFT4»r,AMMA,«t;,yr,,zr. .hfmafx.hfmifx.rmtp.rmpp* 

. iXT.*lP,MPn|NT,MAXPTS,RFLLP,PLNPMm),Pl.NnRi«M?l 

CnMMn»i/RLnCKG/YCRI.«?CRLtHFMARL.MFMIRl. 

COMMOM/MATH/PI ,PIZ,Pin2,nTOR.RTnn 
CnMMnN/PniNTS/NFOGP»MlNTR 
CnMMON/SlIRFIT/NFPOFfT. SIGMA, ftISFAC 
RFAI. MMAT 

INTFGFR SFnGF.SllRFAr.F 

niMPNStON P(5,NTX,NPX),POin(6),PMPW(S),orNT( 5),PBLX<S»,A<^,!»), 

. Rl ^,7 J ,RRn,?J ,MHATm ,C ( ^ , SR ( ^ , F H 3 ) , FR ( ^ ) 

niMFNSinN FXTRAYl l«J| ,FXTRAZ( 19),FXTRAX| 19, 19 ) , FXTR A (1 OR3 I ♦ TFMPH7 » 
OIMFNSinN ZXK 191 ,ZXM( 19» ,ZYI(19) ,ZYM(lq) 

IFIMFFflFlT.LF.O) GO TO 100 
IFIT=19 

on <50 1 = 1, 1 FIT 

FXTRAZI 1 )=FXTRAY( 1 ) =-10 .0*nISFAr.+ 1*01 SFAC 
50 CONTINUF 

Ort 90 I =1,1 FIT 
on 90 J=1,1FIT 
GO TO (f>l .R2,R3,8A) SURFACF 
81 FXT1AX(I,J)= PLNPMTIl) *PLA'nRM(l) 

. -IFXTRAYI n-PI.NPNTI?»)*PI.MnRM(2) 

. -( FXTRAZ ( J )-Pt,N“MTI 3) )*Pl.MnRM| 3) 

FXTRAXn ,J) = FXTRAX( I,J)/PLNnRMl > 

GO TO 90 

87 FXTRAXI I ,J)=-AnRnRF*SORT( 1 ,-|PXTRAY( I ) **?-FXTR A 7 ( J )**7 ) /RFLLP**2 ) 
GO TO 90 • 

83 EXTRAXn ,J)=-SORTt AnRnRF**2-FXfRAV( I) **2-FXTRAZ |J ) **2 ) 

GO TO 90 

84 FXTRAXI 1 .J) = (-4.n*AnRnRF**2 + Fxr«AY| 1 ) **2 + F X TR A Z (J ) **2 ) / 1 4 ,* AORORF ) 

90 CONTI NtiF 

CALL SURF 1 I TFIT.IFI T.FXTRAY.FXTRA/.FXTRAX, IFIT.ZXl ,7XM,7Y1,ZYN, 

. ZXYl l.ZXYMl ,ZXY 1 m,7XYMN,755, FXTRA,TFMP,SIGMA, IFRR) 

100 ALPOAR=ALPHA*nTOR 
KFTAR=RFTA*nTOR 
GAMMAR=GAMMA*nTOR 

A( 1 , 1 )=cnS( ALOMAR )*CnS(GAMMAR ) - S I M ( Al. PH AP ) *S I M I RFTAR ) *S INI GAMMAR ) 
A( I ,?)=SIN( ALPHAR)<-.CnSir,AMMA«)+CnS( AI.PHAR )=!>SIN(HFTAR)*SINCGAMMAR) 

A ( 1 , 3 ) =-Cf'S ( BP T AR ) »S I Ml GAMMAR I 
A( 2,1 )=-SIMI Ai.PHARlsCnSIfiFTARI 
A( ?,? 1 = CnSI ALPHAR)«CnSI RFTAR ) 

A|?,3)= SINIRFTAR) 

A( 3, 1 ) =cns( ALPHAR )*Sl WIGAMMAI l+STMl AI.PHAR ) s-S IM I BFTAP ) *COS I GAMMAR ) 
A( 3,2)=S1N( ALP'^AR)*SI^'IGAMMAR)-CnSI ALPHAP ) *S 1 N( BFTAR ) *COS I GAMMAR ) 
A( 3, 3 )= CnS I HFTAR )*CnS I GAMMAR I 
1F( APROTA.GT.O.O) WRl IF I 20,110) 

no FORMAT! IHl / , lOH***** **»**. 7 X.*ThfTA*,9X,«Y*,9X,*7*,7X,*FRY*,7X, 

. <>FRZ*,SX,*PHASF*,9X,*R*) 


Figure A -2 
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120 


121 

122 

12<» 

126 

12B 

129 


HMTFST-I.OP+60 
NINTR»NI;0GE«0 

PHLKl ^)*P«L<(A)«Pt^LKl5)*0.0 
nn sooo ip^i.np 

OPGPH?tiP|6.1,IP)«!UOO 
IP(APm)TA.GT.O.O) W?ITE<20,120I OFGPHl 
FnRMATdX.opHlso.FlO.Al 
on 4000 IT«1,NT 
t)PGTHPT«P(4,lT,IP)*KTnn 
SINP = SIN(P( S.fT, i»n 
C0SP=CnS(P(5, IT.FP1 ) 

SINT = SlN(P(4,lT,IPn 
cnST=cns(P<4,iTiiP) ) 

B«( 1, 1 )»SINT*C0SP 
R«(2.1)*SINT*SINP 
HHn.ijscnsT 
BBn,2)»-*-FPED(n 
BB<2,2)=+FFEn(2) 

BHO,2)=+FEEn(3) 
call '''ULT32»8,A,B8) 

GO TO (121, 122. 124, 126), SURFACE 
ARaO.O 

BR=rt< 1,1 )«PLNnRM( 1 )+fl(2;i >*PLWn«M12)+B( 3, 1 ) fPLMDRMl 3 ) 
CR=-(B{1 ,2)+PI.NPNT( 1 ) )««PLNn^M( 1 » 

. -(B(2,2)+P'.NPMT(2) )»PLNORM(Z) 

. -(B(3,2)+PLNPNT(3))*PLNn<M{3» 

GO TO 128 




‘ ■'rv 
o 

'ti- * 
r 


Ail = H( 1,1 )*#2/An20RFO#2+(R(2,l )»*2+B( 3, 1 ) ««*2 ) /BFLLP**2 

HRs-Z.OOlrtl 1, 1 l^BC 1 ,2)/AnRnRF«02-Mft(2,l )>!‘R<2,2)+rt(3,l)*Bn,2))/ 

, rtFLLP4’«2) 

CR = H( 1, 2)442/ AnR0RF4'»2+(B( 2,2) *42+8 ( 3,2 )**2)/8ELLP«-*2-1.0 
GO TO 12B 


AR=H( 1,1 )4fl( 1,1) +B( 2,1 )4R 12,1 )+R 13, 1)48(3,1) 

BR=-2.4(R(1,1)4H( 1,2) +H1 2.1)48(2, 2) +R1 3,1) 4R1 3,2)) 

C«=H( 1 ,2)4B( l,2)+rt( 2,2)»H(2,2)+B( 3,2)*H(3,2;-AnRnRF4A0RnRF 
GO TO 128 

AR = H(2,1 ) 48(2,1) +Bn,l)4B (3,1) 

8R=-2.04(R(2,1 )4B12,2)+R(3,l)4h{>,2)^2.0*An«nRF4H(l ,1 ) ) 

CR = B( 2 ,2)4ft( 2.2)+B( 3,2)4 8( 3,2)+4.04Anw0RF4H( I , 2 ) -4 .04A0RnRF**2 
)F (AR.LT.l.OF-10) R=-CR/8R 
IF (AR.LT.I.OF-10) go to 129 
R=(-BR + S0RT(BR4H3-<,.4AR4CR))/( AR-fAR) 

CONTINUE 

XO=H( 1 ,1 )4R-B( 1,2) 

Y0*H(2,l)4R-fl(2,2) 

Z0=fl(3,l )*R-fl(3,2) 
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i ; 

i 

f • 
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IF(NF60FIT.LP.O» GO TO nO 

XOiS««Fn?(YO,ZOf¥MnRM,iNnRM,IFlTt !F1 f f FXTRAY,FXTRAZ,F¥T«*X,IF!T 
. ,FXTRA«SIGMA) 

Rs:SORr( (XO+RI 1,2) )»«-?+»Yn+H(2,?) ) #<'2+ ( 70+R ( ^ ,2 ) ) **2 ) 

XMAGaSORTU .0+YNORM»*2+ZWnRM**2 ) 

NHAT(l)3l.O/XMAG 
NHAT«2)»-S1GN{ YNORM,YO)/XMAG 
NHAT(3)=-S1GN( ZN(WM ,ZO ) /XMAG 
Gn TO 

no GO TO (lAl,n2,13A,13#i),SDRFACF •*»•••♦• 

13’ NHATU )=PLNORM( 1) 

NHAn,?)=PLMnRM(2) 

NHAT( 3)=P|.wnRM( 3) 

GO TO 13ft 

I 32 NHAl 4 1 )=-XO>»BFI.LP <'*2/SnRTU,0>)“S?>>RFLLP'»>!'A+( YO*v2-,ZO**2 )*AnRnRF**A> 
N.lAT(2)=-YO*4nRnRF>!‘>i-2/SORT(‘x'0«»2«HPLLP>!‘**+( YU*#? + 70**2)*Ailft0RF*>»4> 
MHAT{3)=-Z0*AnRnRF*«2/S0RH XO**2*HFLI,P=i'-!‘4+{ YO'»>2+ZO<!'*2 )*AnRn«F*»*> 

Gcr TO 13ft 

134 NHATU )=-XO/AORnRF 
NHAT(2)=-YO/AORnRF 
NMAT(3)»-ZO/AnRnRF 
GO TO 1 3ft 

136 NHATU )a2.0'»AnRORF/SOOT(4.0>.'‘AnRn»':>!‘*?+YO>!>*2 + 70**2) 

NHAT C2 )= -YO/SORT ( 4,n>SAnRnRP=»*?+V0<‘*?+Z0#*2 ) 

NHAT< 3)= -Z0/SORT(4.n<'AnRnRF*w?+Y0«*?+70**2 ) 

1 3ft SCALAR=2.0»<B( 1 ,l X-NHAH 1 )*R( ?,1 )>>NHAT(2)+R( 3,1 )>!'NHAT( 3) ) 
on 1500 I»l,3 

1500 5RU )=B< 1 , 1 )-SCAI.ARANHAU I ) 

FTI=P(l,IT,tP)/R 
FP1=P(2,1T,IP)/R 
C( 1 )=r,nsTpr.osp*FTi-siMP>kFPi 

r. ( ? jacnsT^siNPAFTi +cns?>i‘FPi 

C(3)=-SINT*FTI 
on 7000 1=1,3 
Fill )=0,0 

on 2000 0=1,3 

2000 FKI ) = FIU )+AU,J)4C(.l) 

SCALAR=2.n*( FI { I )»NHAT( I )+F I { 2 )*nwat( 2 ) + F ! ( 3)tMHAT13) ) 
on 2500 1=1,3 

2500 FRU )=SCALARoNHAT( 1)-FI{ I ) 

Y = YO+(XC-XO)=>SR(2)/5P» 1) 

Z = Z0+«XC-Xn)'»SR(3)/SR( 1) 

n=SORT( {XC-XO)<-(XC-XO) + (Y-^Yn)*( Y-Y0) + ( 7-70)*( 7 - 7 . 0 )) 

PHflSF = PI2'i=(R + n)/Xi.AM+P(3, IT, IP) 

PNFWl I )=P«I.K( 1 )=Y 
PNFW(2)=PBI.'U2)=Z 


Figure A -2 (Continued) 
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PNFW(’I|.FR(2) 

PMFU(4)«FR(3) 

PNFW(S)=PHASF 

TFST= HFMAFX*HFMAFX*t<FMIpx*HFM[PX-HFM4FX«'MFMAFX*(7-2C')*(ZrZCI 
. -HFMIFX*HFMIFX»(V-VC)*{Y-YCI 

TF$TBI.=HFMABL*'HFMABt*HFMIBL*HF*»IB«.-HFMARL*HFMARL*(7-7CBI. )*(7-ZCBL». 

. -HFMIBl.4HFMI«L<'(Y-YCBU*fY>YCBL) 

IP ITFST) 2701, 2501, 2A01 

2501 NFDGF=MFDGE+1 FOGF 

SPnGF = MAXPTS-NFr)GF 
IF (TFSTRL.LT.O.O) GO TO 2510 
.U- MnvFM(PBi.K,P( l,SFnGFl,5) 

1F( APROTA.GT.0.0) HRITFr20,2505) PBLK TAPF 20 

2505 F0RMATIlX,«S«,23X,2F10.4,2Fiq.7,F10.4,12X,*FnGF PfllNT, BLOCKFO*) 

GO TO 2515 

2510 CAU. MnvFM«PNFW,P( l,SFnr,F>,5')‘ 

IF( APROTA.GT.0.0) W<U TF ( 20,2512) PMFW TAPE 20 

2512 FfAMAT(lX,*S*,Z3X,2F10.4,2Fl0.7,F10.4,12X,*FnGF PfllMT*) 

2515 CONTINUE 

GO TO 2B00 

2501 NINTR=NINTR+1 IMSIOF 

IF ( TFSTRL.lt. 0.0) GO TQ 2510 
call MnVFM(PBLK,Pl 1 ,NINTR),5) 

IFCAPROTA.GT.O.O) I TF ( 20, 2605 ) PFGTHFT , PBLK , R TAPF 20 

2605 FORYATI IX, <=«#,nX,3E 10.4,2^10.7, 2F 10.4, l&X,«'8LnCKFn*) 

GO TO 2615 

2610 CALL MnVFN(PNFW,P(l . nIMTR) ,5) 

IF(APRDTA.GT.O.O) HRITF{20,26I2) OFGTHFT,PNFW,R TAPE 20 

2612 FORMAT(lX,2t4,lTX,3Fin.4,2FlOi7,2Fl0.A) 

2615 CONTINUE 

2701 IF (IT.FO.l) r,n TO 2«00 OUTSIDE 

IF (TFST*TFSTO) 2704,2800,2500 
2704 ZTEST1=ZC-P0L0(2) 

ZTEST2=ZC-FNFW(2) 

IF( ZTFST1=>ZTFST2.LT.0.0) GO TO 2«00 
CALL INTFRP(POLn,PNFW,PlNT) 

NFOGF=NFnGF+l 
SFOGF=MAXPTS-NFOGF 
CALL M0VFM(PINT,PBLK,2» 

Y=P!NTm 

Z=PINT(2) 

TESTRL=HFMABu*HFMAHL-HFMrBL«MFMIBL-HFMAROHFMAHL*( Z-7CRL )*(7-ZCBL ) 

. -HFMtRL-HFMlRL=>( Y-YCRL )*(Y-YCBL) 

IF (TFSTRL,LT.n.O) GO TO 2710 
call M0VFM(PBLK,P{ 1 .SFOGF) ,5) 

IFIAPROTA.GT.O.O) W? I TF ( 20,2705 ) PRLK TAPF 20 
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270S 


2710 

2712 

271S 

2000 


2’iOO 


3000 

4000 

5000 

502*i 


5050 


Fnft*4AT(lX,<kS«,23X,2FlQ,4,2fil0.7.F10.4.10X,*lMTFAPni.ATPn, Rt.nCKAn«) 

GO TO 2715 

CAU. HnwEM(PINT.P( l»SI?nr,F»,5l 

IFIAPRDTA.GT.0.0) W«1TF(20,2712) PINT TAP# fO 

F0«NATUX,*A«,23X,2F10,4,2F10.7.FlO,4,l0X,OlNTFRPni.ATp0*) 

CONTINUE 

CALL MnVFH(PNEM,P0L0,5) 

TFSTO=TFST 

7FST»<nFGTHET-90.0)**2+(nFGPM1-100.0)**2 

IF (TFST-BMT6ST1 2 9 AO, 3000 » 3000 

BMTFSTaTEST 

BMT=nEGTHET 

KNP=OEGPHI 

HMSRX=SRm 

RMS3Y=SR(2» 

RMSR2=SR(3) 

CONTI NUfi 

continue 

CONT I NUE 

PRINT 5025, N1NTR,NF0GF 


FORNAT»<« NUMHPR OF INTERNAL POINTS.,..., 

, * NIIMBFR OF FOGF. POINTS 

CnSRMTP=RMSRZ/SORT(«MSRX4*2+HMSRY*#2+RMSR2**2I 


*15/ 

*15) 


SINrtMPPaRMSRY/SORT I HMSRX**2+H**SRY**2) 

HMTP=RTO0*ACnSI CnSRMTP ) 

HMPPaRTOO<SASlN( SINHMPP) 

PRINT 5050, RMT,RMP,RNSRX,«MSRY,HMSR2,«NTP,RMPP 

FORMAT!* THE FFFO MAX IS THFTAs*F7.2/ 

. * PHI**F7.2/ 

. * THE RFFLECTFO RAY VECTOR IX,Y,Z) *3F7.2 

. /* THE RFFLECTFO beam max is THFTAs*F9.4/ 

. * PHI**F9,4> 

RETURN 

END 


Figure A -2 (Continued) 
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C 

C 

C 

C 

r. 

c 

c 

c 

c 

c 

c 

r. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r. 

c 

(] 

c 

c 

c 

c 


c 

c 

c 


c 

r. 

c 
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* 7xVm 1 ,7XY !H,?YYMW. tSI.PSw.ZP*TPMP. 

* SIf;w4,TPHU 

THIS SUHRiUlTIMP nPTP'5»*lMP5; T«P ‘>A^AMPTFR5; WPCPSSAQY 
Cn**PIITP AM INTP-lPni.ATORY SilRFAr.F PASSiMi; T^RnuAM A RPCT- 
AM!;ilLAR GRTf) HP FOMCTIOVAL VAI^iip*;. tmf SMRPACF nPTPRMJNPn 
CAM BP RFPSFSFNTPn A*; T^P TFMSOP i>Wnni»CT OF SPI.TMPS liwnpp 

TFM<;inN, TMF X- AMO y-parttai. nF<T>i/ATi \/ f«; ARniiMn tmf 

HOdMOARY AMO TMF X-Y-YARTTAI. oi^UX/AT?WPS AT TM*= FOUR 
Cn^MF^S may hF SPI^CIFTFO nMiTTTFn. for ACT'iai. mAPPTmr 
r|F POINTS OMTH TMF SiKPftC*: IT IS ^iF^. F <; <; A R V TO CALL TMF 
FH««CTIf’M SI WO?, 

MW Input — 


N IS THF nmnkfr i«f i;>in lin^s im tmf x-niPFCTinw, i, f. 

> jMFS parai.i.fi. T'i tmf Y-axis (m .f,f. ?i. 

M IS THF MIIMKFR OF GRIO LI'<FS I"' THC Y-O I « Ff.T I ON , I. F. 

I darallfl ti> tmf X-axis {m .of. ?). 

X IS am array of tmf V x-CnoROT matfs OF IMF oRin LIWFS 
IM THF X-niRFCT I OM. T'^^FSF SHniii.n HF STRICTLY TmCRFASImO. 

Y IS AM ARRAY f)F TH^^ Y-C«’nR n I v a s RF jmf rwin I.T^ifS 
IM JUF X-niRFCTT "M. T''FRF SMfiiiLn mf STRTCTLY ivCRFASTMC. 

/ IS am array (IF THF iJi Fliv'CnOMAL VAi.ilFS AT THc r.Rin 

PHIMTS. I. F. ?(I,J) C'imtaImS T«f fhmcTIomAL VAI.Iif AT 
(X(!).Y{J)) FIH T = 1 AMO J = I,...,M. 

1/ fS IMF RilW r»T*<FMStOM uF TH= •MATRIX 7 HSFn IM THF 
r.Ai.LTNi; drooram i r/ .of. ■-). 

7X1 Amu 7XM ARF ARRAYS I’lF THF ’ X-PARTIAI. OFRI \/AT I mFS 
OF fHP FlllMCTlnM AI.OmO T'R*- X(I) A'vO X(M) ORin LT^'FS, 
^FSPFCTI VFl.Y. IHIIS 7.X1(.t) AMH 7yi.,(.I) COMTAIM THF X-wART- 
ML OF R T WATI VFS AT Ihc pnfN.xs (<in,Y(.l)) Awn 
(X(H),V(.n). RFSMFC I H/FLY , FOR .1 = FITHFP OF 

XuFSF paRAmFTFRS "III hF FOMikfi’i ( /\n I* ApPROXI mAT toms 
Siiopi IFO TNTPRMALl Y > IF ISI.PSM Si» M-njCAlFS, 

7.Y1 AMii 7YN ARF arRA^s "F THF V-PArTIAI, Df.< I \/ot I wfS 
i,c xhf fiimCTIiim ALCWti T"o ''(I) amt. y { ^' > i;u>|o Ll'-c-;,. 
RFSPFCTI v/FLY. thus 7yiii) Amo ci'mTAI'' Thf y-pART- 
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c 

c 

c 

c 

c 

c 

c 

n 

c 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

r. 

r. 

c 

c 

r. 

c 

c 

c 

c 

c 

L 

c 

c 

c 

t: 

r. 

c 

r, 

c 

c 

r. 

c 

r. 


I At. }»iphiwaTIVP$ at T«*- pitr^iTS <«({)>%¥( I H . 

APS>f:cTfwPty, to< r * 

t^FSP PA»AMFTiF<^S WTLt rtF ianfiiiPn (Awn p TI^ATIfiwS 
SMi>t>LlFu INTFsiMaij.y j if ISU>>S'«» SH rwnfCA 'FS, 

'XV**!, Z'<YlM« Aiwn AaF X^Y-PAH|| 4 i. 

HFRlVATIVFS OF THF F«l-*CThi^' M T*^p Fnui^ (f;n«WF!ii, 

(xni, Y(n), (x(-»),ynn* ixui.vcnh* awo jxc««).yiv' »» 

■IFS^FCTI VF»,Y, AMY C»F T^‘p ‘■•ILL » 4 F piMO^Fn 

*:STI'«!«TI‘imS f mi fi.' oau.y ) |p IS'-FSvj sn ImoKATFS* 


cimtapms a s>JITC‘^ : M iMf; wMir.M Hni»MnA>*Y 
nF^IWATlWF IKiFlp.-.AT I Is (KP- -St T FO AMO 
Sttmi.O Ht- FSTI*««UFI' HY r<fS S"H?nMnvF, Tf'i I'*FTP*»*'' TMF 
LPT 

II = O IF /XI IS tlSF-t-SM*i*>L I »-i' |A >'«0 = 1 

I? = 0 IF ;x>- IS I'SF I PI' (AMI' = I 

r-t .- 0 IP /Yl IS IISP!>-SU‘>iJL I PI' (AMl> = 1 

14 = 0 IF /Y'* IS tlSPK-Sll*"ML I PI' CA'-D = 1 

15 = O IF /\Yll IS USP ‘,-SilPFI. IPM 

(AMh = 1 IITP*PI;WTSF J , 
Is = 0 IF VYVl IS l•^^--.-S 1 l^'l*M -o 

lAA'n = 1 f'lwFSwiSF), 
1 / S 0 IF /rYp'i IS IISPS-S'IMPI. I PI' 

I AMO 

1 •! s o I F /XY»"M Is MSI- -.•-sili>i*i I |- 1 > 


OTHFi<w!S 6 l. 

SF i « 

flTHFWMlSF ) « 
nfMPUMlSF I t 


1 i-iTmfkwTSF > , 


JA'MI' = I iiTmfkkTSF). 

IMPM I S',‘-'S"I = n + ?:' 5 ( y + 4 «M + S«I 4 + I A::! IS + •X?#TA 

+ SA* I * + I ;>Hit I M 

CMMS TS'.FSW = 0 iMITlt'.nTP'^ ALL OF<|'/ATt''P ImFHFvaTIUM IS 
f pri ai'M IS'.PS" = I'TIICATPS Mfl nFPTVATIYF 

I MCfV< r I IIM TS I 'SP\-SIII' 'I. T *^0, ".''lY '/Ai.mf hfT'**PFm TMFSP 

LTMI IS IS VAi. lo. 


y,> IS AM A‘?,lAY MP AT i.PnST i.nr.ATIOMS, 

ll-MP Is am A-<</\Y of at l,■'AST •'> + ^'+M l.nCATlOMS wwir.w IS 
ilSi-i) -IM SC-^ATC'* STiHAl.P, 


Amo 


SKivA COMTAIMS lOF Tt^'-STOM FA-TO^:. THIS ''OI.OF IMOICATFS 
IMF r.'i:^''TMPSS nPs{<PO, It- AhSISIu'A) IS mfAmI.Y 7 F«p 

(p, ,oni) T'-*'^ ; r-siii. T > 'F, sikpac.c is afpsox i matpi.y top 
TFMS fp< PsnnncT "p r.oKir spi.Imps, if aosist^-ma) ts i.afi'f 

fi. so. I TMF I FSII'.T I SMt^PArP IS APPROX I wATP'-Y 
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o n .“5 rt n o .1 r. n n o r> r» 


c IP SfAMA pOHAtS TI^SflR l»i|fmilRTS np 

C dime SPI.TNFS ^PSm.T. A STAWft&fltn VAi.(<6 eo»t S)At«tA IK 

C AJPi<nX!MATFI.V 1, IM 4KSnj.llTP »/AI IIP. 

c - ^ _ 

c OM OUTPUT — 

C 

r. COMTAI^S THF VAI,Up<; OF THF xy-, VY-4 AMf) yXYY-PART I A(_ 

C OPR I VAT I WPS OF TUC SURFACF AT TUr OfVPM WOOFS. 

TFRR r.MMTAIMS AM PR^»>R M.AO, 

= n POP MOiMAl. 

= I I F N I «; t.PSS TUAM 7 OQ >» TS l.c.SS TWAM ?, 

= ? TF TMP x-v AI.UPS 04 V-.\*4LI|P«; arf MOT STRTCT*Y 
Twr.^FASTMr:. 

A ^;n 

f*. M, y. Y, 7, |7, 7X1. 7V1 . 7YM, 7YVIJ, ?VYM1, 

7XYVIM, AMD <;friMA A I? p J-pppH. 

T'-'f's SUM40IITIMP R^^•^p-^CMr. p<; packa-;p Mnntii p<; r,pp7, 

A ■'in SNmr.^H. 



r. 

IMTPr.P-^ M ,M. I 7 , t SI.4SU 

< F At. X ( M ) , V ( M ) . 7 { I 7 , M 1 , 7 y } ( M ) , 7 XV ( M ) , 7V ] ( 'M . 7 VA'( M ) , 

7 XY] ) , 7 X V • 1 , 7X Y1 M . 7 YY“«'i , 7l> ( ,M , 7 ) , TPM4( 1 ) 

•IMI = M- I 
M41 - v(+| 

N'M] = M - 1 
MO) r M+] 

MOM = iM+'.l 

1PQ4 = 0 

( M .1. F. 1 -11^. I" .l.c, 1) T'1 

TP ( V ( M ) . l.P . V ( 1 •) 1 (;i» TO A / 

<.»r,MAY ^ AH<; ( S r -'-MA )*P|,MAT { M- M / ( V 1 M )-Y( 1 ) ) 

IP ( ( T .vn. I T <;i.4S'-'/A ) ) ^.n rn ? 

no 1 I = 1 , 

1 7 .J( T .1 ,n = ^Y 1 ( T 1 

TM s 

7 nm.vi = Y( ?)-vn ) 
i)Fi.v7 = n'u.Yl -•■ncf.Yi 
IF (M .r-T. ?) HPI.V? = Y( 7 )-V (11 

IF (HPIVI .IP, n, .- 10 . npi.Yp .IP. OPI.YI) on to a.? 
r.Ai.i. r,pP 7 ( nPi. Yl ,opi,Y 7 , SI o A A Y .r, 1 ,r? ,r.-i .M ) 


tq 




Figure A - 3 ( Continued ) 



nn ^1 » I ,»» 

% « ci*7n .1 

IP {ft» ,F 0 * ?» (;n Tn *» . . 

OM 4 I * i , M 

4 7.p(T,i.n * ?p(i.\,n+f.'^«7(i*^) 

H IF ( n StPS'-'/l*' .''P. ( I r,n TO 7 

nn 4 T » i,« 

MB I s w+l 

t, TFM0(MPTJ * 7Y-»'n ) 

r,q TO in 

7 OFI.VM = Y^^•)-V|M:^l » 
nP|,V^iM = OFI.Y^'+OFI.Y-'t 
IF (W ,nT. ?I OFLYM'* = Y(MI-Vl^t-?1 

IF {OFI YN .I.F, n, .OB, YAt'- ,l.c, nPl.YMl OO TO 47 
Ci'U- CFF7 (-r>F'.VM,-nF|.VMM,^Tf;MAY.C1 

DO Q I = l.M 
MPI = N+I 

» jCMyjMon = Cl*7( I ,M)+C7«'7( I ,MM1 ) 

IF (N ,FO. ?) r,n TO in 
on c> T s j ,M 
MMT = N+T 

o Tcw»(,\iDn = Tcmo( i«io 1 ) ( I , M-? ) 
in IF (yc^M ,i.F. x(l)1 nn Tf» /w 

SinMAV = ABS t S f OMA)-:=Fl.n«T( ) /( Y( '< )-y { 1 ) ) 

IF ( ( I St.PSW/? 1-7 .M*:. '-n TO \? 

nn u .1 = 1 »o 
11 7P( 1 ..I ,? ) = 7X1 ( ,1 ; 

1? IF ( ( T SLP<;'''/^^ )*? ( I ‘^1. 1 ^ 

« ( I <;) p<;w/l )*^ .PO, ( 1 ^;i uSi-i /A4 ) ) no TO 1^ 

OFl.x 1 = X ( 7 )-X { 1 ) 

OFI.X7 = OFI.X 1 +''>*0,X 1 

IF (M .r,T. 7) oFt.v? = x(^)-x(n 

IF (I'FI.XI .I.F. n. .OP. :u-I.X7 J.c. OC| yM On TO 47 
CALI. CFF7. ( nci_x 1 , OF'.X 7 , s in'I'v X . r 1 , c 7 ,CB, 1 
IF ( ( 1 ';i.P‘^'-'/7 ) *7 , t : \ 1SI.P«>M) no TO IF 

1)0 I B .1 = I ,M 

IB 7P(1..I,7) = ci*7( 1 . n+r.7=:=7(?.,|) 

IF .FO. 7) on TO IS 
on 14 ,1 = 1,M 

14 7.0(1, .(,7) = 7P( 1 ..).?)+r,B*M ^,.1) 

15 IF ( ( T SLPSM/B7 )* 7 ,f'F. (TSLPSw/ia)) on TO 14 

70( 1 .1 ,B) = 7XY1 1 

on TO 1 7 

IS 7P(1,1.b; = C1>S70( 1 , 1 , 1 H-r.7:=''.P( 7, 1 . n 

TF (M ,0T. 7) 7P(l,l,B) = 7 0( 1 ,1 ,B)+CB*7P( B.l . 1 ) 

17 IF ( ( I SLPSW/17 • 1*7 ,’MF. (fSL oSM/(S4n Oo TO 1« 
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^XYIWS > 7XYIW ' ‘ “ 

GO TP I «) ‘ • 

H 2XYIWS a C1 <<TGmm( )4.q;»a1-cmip( M4.;> ) 

IP (M ,GT. ?) 7XYIWS s 7XYl^'S+C'=<»TFMDjtM+5i 
IF ( nSLosw/4)#? ,MP. ( isi.PSw/?) ) r,n Tn ?i - 
nn ?o .1 = i,N 

MPMOJ = MPM4-.I 

TFMPtNpMp.n = ZXM(J) 

?l IP ( nsi.PSw/f./*)^? .P^J, ( TS«.PSw/‘»7) 

<• ( isi,psm/?s 6|#2 ,pti, ns».psw/i?«n Gn Tn ?4 

opi.x-v = x(M),x(MMn 
nm.xMM = nei.x*^'+nF».XM 
IP {f* .GT. ?) nP(.XMM = X( M)“Y{ wi-P) 

IF (nPi.XM ,i.F, 0, .nn, npi_xv^’ .LF, nFt.x«*) go in 47 
CAl.L (>"FZ (-nFi,XM,-nPi.XMM,siG\My ,ci 
IF ( nsi.PS'v/4)*p .pn. I isi.PS'^/?) I Gn Tn ?4 
nn ?? j = l,.v 

MPMOJ = MPM+J 

?? JPmP(NPMPJ) = C ( M. J > +C7^.'<7 ( v*MJ , J ) 

IF .FC. ?) GO Tfl 74 
nn 7H ,1 = 1 ,M 

MOMPJ = MOM+J 

7H T = MP( MPMOJ) = TFi*o (MDvpj j ) 

?4 IF ( ( T SLPS'J/ft^* ) *? .'-'F. ( I SI. PS'J/^7 ) ) GO TO ?*; 

= 7YVM1 

GO Tn 74 

?S 7P{M.i,^) = Cl =^7P < I , I )+C?=:=7p ( MW) , ) , 1 ) 

IF (M .GT. ?) /P|-i.l,^) = 7 o| M . 1 , ^ 1+G:<*7P ( M-P , 1 , 1 ) 

7S IF I ( I si.PSw/ps'-> ) *7 .'IF. { I si.os'.'/i ?« ) ) <;n to ?? 
7XYMNiS = 7XYMM 
Gn Tn ?° 

77 7XYN''IS = c I* T F*o I MO-; ) +r. 7*TF MU ( MpM_ 1 ) 

IF ( vi .GT. ?) ’XYMM*; = 7XY'"MS + C'a*TFMP( ) 

nO| 1 = Y( 7 )-Y n I 

IF ( GFI. 1 . 1.F. f). ) Gn TO /.7 

npi, I = l./nFi. 1 

nil T = i»M 

7G 70(1.7.1) = nr-i,T*( 7. ( 1 ,7)-7 ( 1 .1 ) ) 

7.P(1,7.H) = nFl. («(/■>(!.?,? )_7p{ 1 ,1 .?) ) 

7D(Ni, 7,:<) = n=i. r<( TFMP( vUM+>' )-TF'.'P{ MPf,'+l ) ) 

CAi.i. TFi^MS (n» AGl ,SnIAG) .SMv'AY.nciJ. ) 
nl AGT = I. /ni.'.GI 
nn i = 1 

TO 7p(T. 1,1) = nl AGI=M 7P{ I , 7 , 1 )-7p( 1 , 1 , 1 ) ) 

7.o(i,i.-^) s f'f AGI«{ 7P( 1 '^)_7p( ] , ) ) 

7P(<-‘,1,'^) - nlAGI#( 7 P(m.2,'X)-70(m.1 ) 
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1 ) s n|*(;i*sr»? API 
IP (w ,P0, ?) r,n T'» 
nn j » ?,NMi 

,IMl a j-i 

J»l « J4-1 

NPKtOJ a MPM+J 

npi.? = Y(jpn-vu» 

IP C»pl;> *».p. n.) r,n TfJ A? 

HPU = 

no I = i,>* 

^pn,.joi,n = npi.i*f 7fi ,jPu-7( I ,ju 

a riPl.I- (7.P( 1«JP1 ,?)-7P(1 ,J,?n 
7P( “♦JPl ) = nPI.I#nPMP{ •'•PMUJ+l )-TP*'P(MDMPJ) ) 

CAI,I. TP^MS ( ni Ai'?. sni Ar,?,si';MAv,r>P|j) ) 

niAniM = I , / (O! AGl + ni An?-sm A(;i *T‘=‘^P(vlMl ) J 
nn I = I.v 

7P(T,J.ll = ni AGin«{ 7P( I ,.|PI ,n-7PU 

* Sni a ;i *7(>n , , 1 ) ) 

7P{1,.I,-^) a n? AGI ( 7p ( I , |U)) )-7P(l , 

* Snj AG! *7.P(1 ,JM1 ,^) ) 

7P(w, = o! 6f,I 7P( JPl ,•:< )-7P{ M, j 

* Sni AG1 *7 P( -n.lMI , :< ) ) 

TPMP(J) = n| AGf M*SnT AG? 

ofAr;i = niA';y 
si*TAf;i = S'UA(;x 

^4. niAG fM = I , /( nr AG j-sny AGi ❖T'= -p ( M-.(i ) ) 
nn HS I = i,M 
NPT = I'l+T 

HS 7P{t,M,l) = nl T'-'^PI '■PI)-7‘J( T )- 

==' AGl := /p( T . V‘-M , I ) ) 

7P(1.M,H) = OlAGlM*(7XY1MS-7lM 1 
« S"l AG\*7on .M‘^1 ,? U 

|PMO{M) = nl A('.fM*( /XY'''";-7P( 1.^)- 

« Sn I A(, 1 :;:7 P( >1 , ^.M 1 , ? ) I 

nn ^7 ,1 = ?,M 

,mA< = mo^.,1 
.MAKPI a ,)H.A<+1 
T = rP~»P( JMA< ) 
nn A A I = 1.1^ 

7P(T,JHAK,n = 7P ( I ♦ ,)HA< . n-T«7P< I . JMAKPI . n 

I .>MA<, ^ ) = 7«-’n ..IHA< , 7 )-T«7P{ 1 ,.|HAKP1 

H7 T^MP(.IHAK) = 7 P ( M, .|KA< , ^ )-T*Tc ( ,|H A < P] ) 

M'^U = X(?)-X( 1 ) 

IP ( npl.l .i.P, n. ) <;n jn 47 
nn. I = 1. /’■•='. 1 
nn ■•$«? J = l.M 
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?*»(?. j,?) = nFi.T*!‘( ?J ..m 

CAU. TPilMS (nlAr,l,SI'TAai,5iTf;MAX,nFl.n 

niA^,l = i,/niA(;i 

nn j = I ,w 'Uniff , 

7.P(l,J»7) = IM A.;i#« 7P( 7,j .7i-7pn ) 

HO 7.0(1, j, Hi = nTAr,l#(/o(?,J.H)-7P(l, .1,^)1 
TFMKM+ 1 ) = ni Ar,i#sni Af;i 
le (M .FiU 7) i\i\ TO 4H 
on I = 7,0 « I 
1t>^l = I-) 

I01 = 1+1 

NL>I = ^t+( 

OFL,' = X( IPU-X( 1 1 

IF (0 + 1.7 .L+. 0.1 <>0 T(' 4/ 

OFI. 1 = l./l'» + '.7 

ht) 4.0 J = l,M 

7.0(1 FI,. 1.7) = nc|_l«( 7 ( lo] ,.))-/( 1 ,.m ! 

4() ;o(I01,J,H) = 0*-».?«(7Pn^M ..1,1 )-7o( I .0, 1 ) ) 

C.AIJ. TFKOS (olAG7,Soltt(-7,Slt; ’AX,nFl.7) 

0TA(',iM = 1 . /( ni A',;i +ni A<;?-soi A(;i>.-tfmp( MPT-1 ) ) 
nn 4.1 j = 1 ,M 

7.0(1, J,?) = OIASI ''*(7P(T.>1 ,.I,7)-70( T ,.1.7)- 
s< <i''i iM;mt/i->( 1 M) ,.1,7) ) 

4.1 70(1,. I, HI = OT a;; fN<s=( ■/ 0( ! O 1 ,.t, H )-7o( 1 ♦.), HI- 
SS SOT m;)s.s/0( IMI ,J,H) ) I 

TFV.J(Mpf) = (T T a;;i \USS.-»T AO? I 

DlAOl = OlA'.'.P 
4? sniu(;i = s>>l <\C.? 

a-i l)lA(;i'\l = 1 . / ( I'T (\i;T -' ifti A01 =STt-"i>( mu.'-T ) ) I 

nn 44. ,) r 1 ,M I 

MPmP.I = MO\s+.l 

7P(''1,.I,7) = 'M ai;i .:=;t 1 T'- '0( Ii>, ))-/.>(.',.),? )- 

NITT A<;T sw ( V] ,7 ) ) 

44. 7 J ( V,, ,1 , ■:< 5 = ‘M A,;i -.s;t ( 1 .1 ) _ 4 .s ( M, ,1 , _ 

--S S'M AO, 1 */P ( VMI . J .7 ) ) 

f(i| 4S I = 7,^' 

IHA< = noi-f 

IHnK'Ol = THAS + 1 1 

\iu ] >■( AK = '1+ T hA< 

T = rF'>'0 ( MO f AK ( 

nii 4.4 j = ] ,I\, 

70 ( 1 hAK ,.1, 7 ) = { 0 ( I OAK,.!, 7 ) -T*70( MAKP:. , J .7 ) 

4h 70( 1 HT A< , s ) = /0( (HAS,.),H)-T«7 o( TWAKPI ,J.H) 

■<oTn< .1 

44 IFPP - 1 
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piincttom SliRFn? (XX,YY.ZX«7y,M,M,y,y,7*!7*7P^Sir,wA) 

c 

C TMIS PMNCTinN IWTP-^PnLATPS A S'UPACP AT- A RIWPN CnnRniNATP 
C PAIR USTMR A rtl-SPl.U'P UNnPR TPNSI»tM. T^P SUHROiITIMP S'lPPl 
C SHHiltn HP CALLPn PARUPR TO npfp\«lWP CFRTATM mjiCPSSARY 
C PAHAMFTFHS. ’•* 

C 

c nv IMOIIT — . , ** 

C 

C XX AK‘1) YY CflNTAlN THF X- AMH Y-Cnnwni MATFS np THP PHIMT 
r, Tn HP .VIAPPPO niMTU TMP I^•TPRPMl.flTI^''i SMRPAf.P. 

C. 

C. H AMO M CONTAIM THP MIIMHPK HP iV< 1 0 I. IMPfi TM TMP X- AMIt 
C Y-nIRPr.TI ONS . RPSPPCTTx/pi.Y, (IP THP PPCTAmRHLAW r,Rin 
C HUTCH SHPCIPTPil THP Sn<PAr.F. 

C 

C X AMI) Y A»P ARRAYS CHMT A I M i Mf; THP X- AmH Y-fiRTO WAI.MPS, 

C RPSPPCTT YPLY , PACH H' ImCRPASTmc; n«ncR, 

C 

C 7 IS A MATRIX CnMTAlMiMG Thp m =:> m PiimCTTHMA). VAI.upS 
C cnR^FSPnMniMG Tn TMP si'in waijipp (T. p, 7(T,vI) IS Thp 

C SlIkPAC.P WAIJIP AT THP Pfil-MT (X(!).YU)i PPR 1 = 1 m 

C AMO J = 

C 

C 17 CMMTAIM.S thp RHW nlMPMSlOM np thp ARRAY 7 AS nPCtARCO 
C 1 M THP c A l_ I. 1 M(; P R ()(', i A . 

C 

C 7.P IS AM ARRAY OP R*'Vi*M I.OCATTOmS ST.nRPO WTTH THP 
r VARIOUS SIIRPAC= OPRIVATI'^P Impor^ATIOm dptfR^IMPH hy 
C. SORPl. 

C 

C Amh 

C 

C sigma CONTAIMS TOP tpmsIom PACTOR I ITS SIPM IS IPMnwPn). 
C 

c TOP PARAMPTPRS M, M. X. Y, 7. 17, 7P, AMO SIPMA SHOOI.H HP 

c TM-HIT ilMAi.TPRPn PROM TH= noTPoT OP SORPl. 

C 

C OM ooTPOT — 

r. 

C kO'RPo? Comtaims thp I mtp RPoi.at-o <;m'RPACp vaijip, 

C 7X IS PARTIAI. OPRIVATl\/P ''ijo PSMi^CT TO X. 

c 7Y IS pari I At. nP<TVATI\/P HlTH RCSPPCT TO Y, 

C 

C MOMP OP THP IMoilT uaRAmptpRS AR p AI.TPRPH, 

C 

Q IHIS POMCTIOM RPPI-RPMCPS DACKAPP* '.VlOOl.PS TmTRV. Amo 
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n i-> i*» 




HF»Mi) ( FI ,Foy) s 




. allM «P •»>•••> ao^o •>• •aa>v»>aa>x*»»—«>oaBi—^» o^aiXSMilMiM^xM^kCkttAitoabakaA Ik 

fMTFr,F!» M,N,T7 

■^FAI. yX.VYt7.X.?V,ytM|,Y(M),?.( T7.M) 

HFWMZ ( FI ,F?,FDl ,FP?) * ( F^PPFU +Fl /n<F|.H*»nFl. 1* 

hF(.7>i«t ca?«(ncij+nFI.^) + 
poi#(np|. F+PF1.S) )/ 
(A.-XPPI.S» 

(P?-Fl 

■f I ( 7.#nFi.i*nFi. i-nPi.7>!i{ iiFf,.i+nBi.«; ) )*fp7 
-( ?,*nFi.?»tr>p|.?-nn.l #( obi. ?+ n«:i,s ) J^FpI ) 

/f(^.*nFLS) 

( FI .«=7,FPl *FP7.STr,>iAO) = (P7*rtFI.1+Fl>!thBI.7)/PI!».S 
+ (Fi>?*( <;TMMMi*nci.?-nPi. l*( 7. * f r.nSWPI + t 4 )♦ 

<; I^MOP + S Tra*^AP*rn^upl>ithPl. 7 ) ) 
+pjl*( <;tMHM7#npi.i-PFi.?«( 7.#(rn«;MP7+i * )« 

<; TMM01 +«s ICAO# f.n«;MP7*nF 1,1 i ) 
)/{ srr,MAP«SlfVAAP>i<nPi.s«( siM>-'f^«: + SK'MAp»itnB(.<; j t 
MPPMK'O (FI ,F?.FD1 ,pD7,mf;MAM) = ( p ?_c n /nPI. <; 

+ ( ( np».s*STG‘j'AP*f:n«;Mv.| -SI )«PD7 

-( nPi.S*ST r,M AP«r.nSM'<?-s I mmms ) «poi ) 

/ ( S rn'KAo#*;; ri^M AOs:!npi,S '‘( S ( + s f r^MAPftopi^ <; ) ) 

= AHS( ST';M^ ) »:.ci nA r ( M-1 ) / ( v( M)-V (1 ) ) 

= AHS( SI r,MA )«Fl.nAT(M-l ) / ( V(^•)-Yn ) ) 

AM 

,>- ) 


SIGMAV 

SIr;^'AV 


,|M1 = 

tMHVI. (YY.V 

,1 = .1 

1 + 1 

IM1 = 

(yy.x 

T = 1 

.11+1 

riFi.i 

= VY-Y( vlA 1 1 

OP).? 

= Y(J)-VY 

nci <; 

= Y ( .1 l-Y ( IMl ) 


|C (s tom AY .A't-. .).) an Tn 1 

7.1M] = ( 7 ( T-1 ..l-n , / ( I -1 ..I ) . 7P( I-l ,.1-1 . 1 ) . 

70 ( 1 - 1 . I. in 

n = MBUM7 ( 7 ( I .1-1) .7 ( 1 ,.l ) , 7 J( 1 . , 1-1 . 1 ) ,7D( I .,1 . 1 ) ) 
7XXT''1 = MFPm7( 7P( 1-1 ,.|-1 ,7) . 7i)( T-1 .1.7) . 

: 7P ( T-1 .,1-K -1) , 7P ( 1- 1 ,.| . 3 ) ) 


7XX1 = ( 7 II ( I ., 1- 1 , 7 ) , 7 U( r 


7 ) 


1 ).7P( T-1,.1,1 ) ) 
j. n ) 


- 7P( I .,1-1 ,7P( I ) 

/YT«'1 = MPpMn ( 7 ( l-I . .1- n , 7 ( r -1 , .11 , 70( I -I , .1-1 

7Y1 = ( / ( T ..I- I ) . ^ ( 1 ..I ) . 7P( T ,.1-1 . 1 ) » 7P( T 

7VXV(v. = MrPMH ( 7p( T -1 , .1-1 . 7 1 . /P( I-l .,1. 7 1 , 7P{ T-1 ,.1-1 , •^ 1 . 

»!< _p(i-i.,i.^n 

7YXX1 - MF-Xiv^n ( 7P( T, J-1 ,7),70( I,j,?),7P( T..I-1 ,7).7D( 1,.1,-n ) 
r,n Tn 7 
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I 0FJ.P1 s (DPLl+nci.s)/?. 

OPI.u? = (l)F(.?+nP».S)/?. 

r.Al.l SNMCSH ( STWHMI ,Sf(i'*AY>«'nPl.l ,01 ^ 

r.Ai.i. SMHCSw (’>IMMM?,cnsHM?,sTr.<-Av*nP(.?,0) 

f.Al.l. SNHr.S^^ ( '»My<;,n>lMMY,ST<4«»&V«np».S,-n 

CAIJ. ( SfMMpl .nilMMY.ST^;'''' 'ftOPt-.l /?, ,.-n- 

CALL ( <;fMHU?,nMMMY,STUMAY’!‘nPL?/?.,-1 ) 

call .SNHCS« { riMMMY.cnsHPi , Sf <Av>!<nPLPl . U 
I ' LL *sMHCS 4 ( nMW»''V .r.fiSMp? , <4 Tf;MAy#nCLO? ♦ 1 ) 

7IMI = HPKMia? (/ n-1 1-1 ) ,7( f-1 ,.i j ,7P( 1-1 ,.1-K n . 

« 7P( 1-1 . j, 1) , <;Tf; 'AV 1 

7 1 = KP7-IW7 (7( f , J-I) , 7( 1 ..I ). <P( I , J-1 , ^ 1 . 7P( T,«J, 1 ) , 

* '^If;vAY^ 

y.xxiivil = HpyM,.i7(7J(r-i..i-i,:j),/pn-l,.),'>). 

* 70 j 1-1 1-1 .7) ,7P( 1-1 ,sir;MAV) 

7xxr = MP^MM7( 7t>( 1 ,.|-1 ,? ) , 70{ T ,.l, 7 ) . 

* 70{ 1 ,J-1 , 7u( 1 1 ,s ir,MAY ) 

7YIM1 = 7(1-1 .J-1 ) ,7{ I -1 ,.n ,70(1-1, J-l . n . 70(1-1 I, 1 ) . 

* SICMAY) 

7Y1 = HPRMMni 7 ( 1 ,J-n ,7 (1 ,.n , 7 D( 1 ,.1-1 , n , 7P( I ,J,1 ) . .STr,MAV) 
7YXX1 M = hprmv,i( 7i>( 1-1 .1-1 ,70( 1-1 ,.l ,7) ,7P( 1-1 ,.1-1 ,R ) . 

70 ( T-1 ..i,:^1 AVI 

7YXX1 = HPI-'W.1( 70{ r ,.|_1 , ? 1 , / o( ( .,1,-7 1 , 70( I ..1-1 . 70( 1 . ) . 

=:> SI (4 MAY) 

7 1)0). 1 = VX-X( I''| ) 

1)01,7 = X ( 1 )-V < 

l)PI> = X( I )-X{ I'll 1 

IP (SraMAX ,MP. n.) ''.M Ti: 

<4l|i,'P()7 = HPRM7 ( 7 I n .7 I , 7YV 1 "1 , 7 vy 1 1 

/X = MPI^-II) (7rvi,/r./yyrMi,.'fxl) 

7Y = -^1PPM7 ( 7'' T VI , 7 VI , 7VYY 1 ^ . /Yxy t ) 

^ P 1 1 1-.' 

■*. npl.Pl = ( nPl.l + nOi.s 1 /?. 

0PL=^7 = ( HPL7 + .1PLS ) 77. 

CALL SMMCSM ( «; Imwm 1 ,r.nsp;vii , mf;..-AX *nP|, 1 ,0 1 
CALL S'-* ( S T MO.-A7 , Cf.'^M vip. ■; r r; 1 A X=.:0PL7 . 0 ) 

CALL I MM.K ,nif.ivY. S Ti;-.'AXs7nPLS »-n 

CALL SMHC.SH ( S) vMPi .mi.MMY, sr.;.MAY!:<MPi.1 77. ,~1 ) 

CALL SMHC.SH ( SJMH07 .niiv.Mv,si!;vAx*nPL7/7, ,-1 1 
CALL SMHC.SH { ntivMY.cnSHOi . sr.:'-.Ax><ncLP1 , 1 1 
CAI.I. S^‘HCSH ( M'l VYY .rnSMO? , s 1 i-MAX ftDPL o 7 , 1) 

SMRPIV s HPRMM/ ( 7 1 Ml . 7 » , 7 XXT 1 . 7 XXI . s ICHAX ) 

7X = MP^MMfi ( 7 t VI . 7 1 . 7¥X IMI , 7XX 1 , S ir,vAX ) 

7V = Mpj\iM7 ( / V 1 Ml . 7 Yi , 7YV yl ,1 . 7 VXX I , S ICviAX ) 

p P TMO M 

PMO 
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nooof~' rior'tOoOkTOoooo 




SiwBfMJTlNP SNMC^M ( , 1 SW) 

THII^ 'SUM*<'»M>TtN^ s^I^TJKMS APPRHK I .^AT TH 
StNMmX» * SlNH<X»-»< 
r,r>SHiA|o • cnsHiXJ-l 

AMD ■<! 

cns«MM(x) ■ cnsHi 

UITH APLATIVP ».PSS T^AM X.A?P-l<k 


Iiw INPUT 


X C^MTAl^>S THP VAIJIF IIP TMP I'Mnv-PPMnPMI" VAR1AW|.P. 


ISW iMiiiCArPS TWF pumCITum nPSl'^pn 

- IF iiwLY SlNH'^ IS UK<;i.<Fn, 

a 0 IP rtUTU SINH*'^ ANi> CilSHM APF PPSI«PU» 

r 1 IF UMI.Y CnSHf*' IS ''»FSlKFn, 

- 2 IF tiNI.Y i'USUM»‘; IS M^S^*<F’''♦ 

= I F ►<{)TH Sl''‘’H'v^ A'"!’' Ci'SMMM APP OFSIRPH. 

r, 

c ntvi nuTPUT — 

r S1i\--»'*I r/F"TAlNS l^P VAI.UP np SU'Pf'tX) IF 1.SW *'-F. 0 

c ISM :i (SIWMm is llM.vl.TF•<Fl^ IP ISW .pn.l UR TSW 

V. 2 ) • 


HP 

,pn. 


C 

C 

C 

c 

r. 


CnSM.t CUmTAIVS Tpp YAI.IIP ' tP CmSh-U “ ) IF 1 Sw 
I’sw .F-. I •I'FvtTAl «'S T 'P \Fit.uP UP CnSMMN.( 

,\')P, ? (CF'SHN IS ni'‘AL IF ISW .FU, -1). 


,FO. 0 HR 
X) TP IS'-' 


i: /v-'-t' 

X ,\'U> ISW AliP 

c 



!MT.'SP( ISW 
PiU. S I NU«i .CHS'- «K 
UAl A S.->4/^,S02V<’.S'^^S • 1 S'A'^F-W ' / , 

S W 7 S 7 W}‘' 0 ?C- >.u/, 

StM I>U(1F-1'2/ , 

•K S'.'l /-P. SA 4 U' I 7 s I I(IP- 0 '-< / 

OAT t' r.PA/ 1 . 7 »A I nS 474 Wni w(',P--i 7 / , 
r,RS/ 2 ,S 72 ^ 72 P 77 '-J'U) 4 /.F- 'S/, 

r.p;/?, isi sis\‘)Fo?u 2 '-‘"-'>a/. 

S-- C:' I /7 , 1 ■' 1 ^ ' S-'/s4P-'/X / , 
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OATA 

tt 

XX S 
dX = 

.XS = 

Ii= MAX ,GF, ?, 70 ) ,nt?, (flx i.is .awn, 

* rSM ,WF. ;^)) PXPX = f-X-»(aX) 

If' nsw ,FO. I ,(U, ISV) .pn, ?) j;r. to ? 

IF (a;< .gf, i.iG) TO ) 

SINHM = ( ( ( ( ( SP4-XS+SPM«yS + 5;'J?I#XS+SPl )«XS+1. )*XS#XX) 
=!= /( (SGl*XS+l. »VA. ) 

GO TO 2 

I SINHf. = -( { ( I ./PXPX+4X ) + AX I--XPX ) /?, 

IF <xx ,LT. O.) SM'Hm = 

'/ IF (IS-'J .-'tF. I) I^'J .v-. n f;r» yn 4 

IF (AX .GF. l.IG) >;'• 7'> " 

cnsM". = ( ( ( ( (t:i-> 4 »xs+r.;’'o^.= xs+::p? )*xs+CPi >«xs+i. ix'xs) 

* / ( (Cni*xs+i . ) 

G'i Til 4 

T (jiiSFiM = ( ( 1 ,/i-XPX-P. )+FXPX) /?. 

4 IF (ISt-i ,LF. 1) 

IF (ax ,GF, P./n) Gfi TO s 

CnSHi«. = ( ( ( (7“>XS-XS+'Pp)5XXS+'iJ? )*XS + I. )«XS*X*; ) /( ( ( ( I/f)4 

*XS + Z >'H>=?Xs+ 70 >' )*xs + /<il )*xs+l , 

^CT(|^,V, 

5 C(ISH*« = ( ( ( I . /FXPX-7. )-X<; ) + PXO> ) /7, 

FiMP 

S(iHUiMTli\'F tf<-is (nl ft-, sola sigma, mfi.) 

C 

C TMIS SliM-('iin IMF CMwpIlTFS THF nlAG''’''Al. Af'O SHPFPOIAGnwAL 
C TF-TMS ;iF IHF T‘<tnT,}r,f|tMai. | j mc/.^ SYSTFm assnciATFO WITH 
C SPLIiMF UMOFi? TFMSl'l'i iMTFPP'H.ATf 
<: 

c flM (MPIII 

c 

c snThia coNiaiMS thp tf'''smm ‘'A'Mn-', 

c 

C AMI) 

r. 


2?ii/5.‘i<»?'»YU#>?44/?0F-n7/, 
?P?/l,77T4^4‘i‘'»iT^{T«tt4F-n4/, 
7.PI/1 ,A'JAO04(SIF<»47 <)?f-»>?/, 
204/1. ■ 
1 0 ? /*‘5, *T«4^»^'744 1 4 IstM-oy / ♦ 

. 27H 1 494440^4A,hF-=> 4/ , 
ZOl T7 14^ 51 «IF-I.?/ 

X 

AmS(XX) 

xx#xx 


mfmm (imm 



•r 
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n n ri r-, r> o n n n o n n n " o o n n r» r> o o n o o n n r» r> o r> v-> 


J>P«- unNt»lMS tMP StPP 
iiM nUTPMT— ' 


n 


(SlfiMA*npi.4tcriSH(sfKM&*nci,» - siwM(smMA#nptl 

( SlU'^A>S'OPl. * R?MH|S!f:MA#nPI. ) 

SIGMA#nP|. J - STGMA«nPL 


sni&r, = nPL*!'- 


{ ^TMW{ STG«&*hPl.i 


AMO 

SIGMA Awn IIPI, HMAl^TPt^cn, 

TMIS SUHOnitTlMP "Itpc-^CMr.FS PACKAPP Mnmii.p S^^MCSW, 


t?PAI. ni AG.Sni AG.SIG'*A,nci. 

IF {STGM& .IMP. n.) r,n TO \ 

01 AG = OFI. /-^. 

SOI Af; = nFi,/6. 

^FT'R'M 

I SIGOFI. = SKiMA^SOPt. 

CA«J. GNiHCSO ( SI MHM,roSMM,Sl';nP|.*0 ) 

DPMOM = OPl./( { SIMH»<+S100t-I, >*SIO^PL«SIGnPI. ) 

01 Ai'. = 0FMnM<!( siGrici.#r.nsuM-.si ►^p^. ) • 

SniAf. = DFMfVU-tSlMMM 

•<PTII><M 

FMO 

SOHsnoTlMP CPP' (0=1,1 ,npi.7,sn;viA ,i:i .rp.ct^.M) 

THIS snB'<noTlMF OPTP'IMIMPS THt- si'tPFTr f pmt'; Cl. C?, Amo 
OSPO TO OFTF^MImF umopoimT SL0P''-S. SPPCIPir.Al.LV, TF 
P'tMr.TIOM VALOFS YI. Y>>, AMO A<P GiVfPM AT UOTMTS XI. X?. 
ANHV X^. =!FSPFCT I VFl,Y , TOP OnAWTlTV ri»:'YI + r.?*Y? + f.'^^'Y^ 

IS IHP VCAI.IIF OF THF np^^fWATIVP AT VI OF A SPl.lMF IIK-OPR 
TFwsinM (WITH TFMSlMM FA^.T^^ Sir„iAl PASS (Mr; thwOUGm Tmp 
TH< = F POIaTS AMO Mio/IMG thMO o=<I\/ATI''P FOIIAI. TO 7F«n at 
XI, OPT I OMAI.I.Y , r)Mi,Y THii ''AIJips, Cl Aoo C? ARF nPlFPMTMcn, 

OM IMPOT 

OFU I S XP-XI ( ,GT. n. I , 


Figure A - 3 ( Continued) 



n o o n ri n v"> o n 


c 0F«.> IS x^-xi (.r.x If M ,po. ?, TMis PkPAttpim rs 

c laMrmfo, 

c 

SIfi»*A IS THP TPNSin*a f ACTOR. 

Awp 


w IS A SVTTCM IwnlCATlMG THF MiiMHfi# nf COFFFICTFWTS TO 
M*= ^FTllR^FO, TF Ai ,FO, P OMtY Two r.nFPPTC T fmTS ARF 
RFTtlPWFH. HTHFpujSF AM. THPFF A«p t^FTMP^^Fn, 


HiM f'llTOHT — 


C 

C 

C 

C 

C 

C 

C 


Cl, c?, ANn THF COFFcICTCMTS. 

vnWF MF THF IMOHT PARA?*FTFPS A-i F Al.TFPFn, 

THIS SIIKRiuITTNF •<FFc*<CNr_p<; PAf,nA.;c mo’^im.F S^HCSH. 


RCAi. DPI. I ,PFi.?,sTCr4A,r.i -.c?, 

IF |M .Fp. p) r,n ro 7 
IF (Sir, A A ,NF. o.) ',n j'\ ] 

HF|. = |)F|. P-OFl. I 

Cl = -( OFI. I +PFI. ? ) /( I’tFI. 1 *nF!.7 ) 

C.P = nF|.?/( PFI. i#nc( ) 

= -DFI. l/(nct.7«f't:i. I 
• F T I H M 

CAM, SMMCSH {niiMMY,f:nSHM| ,SFf;'Aft<=nF|.i , 1 ) 
SHHCSH { nii .^v,Y .CnSH -M , S T 7 

= STr,N'A*n)Fl.7+nF|.l ) /7. 

= SK',MA«( r>c:'.?-nFU I /7. 

SMHCSH I STM‘-»,v»P.Pt|>.,.iY,nFi.P.-l ) 

S MHC SH (SI , lit |..,r-Y , O F I. f - 1 ) 


CAI.I. 

DFI.P 

OFI.A 1 

CAL'. 

CALL 






1 ) 


OFMn.,1 = cnSHMi *( fipi ?_npi_ ] ) j | s i „iMMP+nFLH )* 


( S!MM^t« + ,■ 1 P| M) 

Cl = 7, *( SINM^’LJ+nFL •*)*<'' f''‘H'' '+‘')FL''^ ) /PFMriM 

C? = -CnSHM 7 /f>F'>iini 


CT = CnSHM] / hfm'Im 
KFT'iPAr , 

7 Cl = -K/OFI.l 
C? = -Cl 
IFTIKM 
FMO 


C 


ciiwCTinK' IMTP'H. (T.'^.M) 
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o n o r> r> o o 


TH15i PUK-CTIf»N nPTP'<'^TM«^S THP np ywp JmTP«VA|. - 

HY a niYPN' Twr<PASlMt: SPOM?:a<CP) JM WMICW 
A MVPN vAi.iip I.IP*;. 

HM IlviV>HT 

T IS T"F aiVPM VAI.IIC, 

C X IS A v'Pf.Tm np sr<ir.T«.Y iwC'-'pasim*; ^fA^.nPs. 

c 

r. AMii 

c 

r, M IS T I.PNGTH HP X (M 

r> 

r, nn niiTPi'l — 

r. 

r. ><HTtHNS I'M 1 «iT>-api- i smo imat 

c 

C I = I IP 1 . 1 . p . X ( 7 ) , 

C I = M-l IF Xt v'-n .'.p, T 

C SP X(T) .l.p, T .I.P. X(I-^)). 

r. 

c >v r np TMF input ‘^A'^ a ik-TF'^S AKP AITp-'F’"'. 

c 

(; 

l-\jrF''F< (\! 

-,Fal. i',x(M) 

TT - I 

IF (IT .l,P. X( -'>) ;•> Ti' A 

IF I i r .UP. X ( ■■ -1 ) ) Tn S 

M. = P 

I ^ = .,s- 1 

1 1 = ii.+i PT M P'.‘>A f ( I f-iUF ( I r-v ( !.. n / ( v ( iH )-x ( ti. ) ) ) 

IF n r .LI. X ( n ) ■ i ■ ' p 

IF ( TI .LF. X( t + ! 1 ) {;•> T" ' 

I'. = I+l 

fv'! lit 1 

p T -• = I 

•;tt lit 1 

^ I ■ I I'M. = I 

' r- I Ms \ 

/. I \t I V/ 1. = I 

' F T' K . 

S I t,' f ;\M. = -1-1 

■ PTl’-..' 

F I . 
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